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Abstract 
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1. Introduction 



In ref. [|IJ, the possibility of determining as to high precision from the hadronic 
decay width of the r lepton has been put forward. It has been argued that, in spite 
of the small mass of the r, this determination is perturbative, since in the contest of 
the operator product expansion and QCD sum rules® non-perturbative corrections 
to this process can be estimated to be small. This is essentially a consequence of 
the fact that in the chiral limit the dominant non-perturbative corrections to this 
quantity behave like \jm\. 

This determination presents several subtle points. First of all, it relies on assump- 
tions that are somewhat stronger than the usual assumptions of perturbative QCD. 
One should in fact assume that it is legitimate to add to a truncated perturbative 
expansion, which is essentially an expansion in inverse powers of logarithms of the 
momentum scale of the process, terms that are formally suppressed by powers of the 
momentum scale itself. A justification of this procedure can be found in ref. ||, where 
its validity is related to the position of the singularities of the Borel transform of the 
perturbative expansion. There it is also made clear that the underlying assumption 
of this procedure is that the only Borel singularities which are present are the known 
ones, i.e. infrared renormalons and instanton singularities. 

Recently some authors!! have brought arguments in favour of corrections, due 
to ultraviolet renormalons, suppressed by two powers of the momentum in current- 
current correlators. These corrections would give rise to terms behaving like \jrr? T in 
the r hadronic width, thereby spoiling the analysis of ref. [|TJ. 

In the present work we will not deal with the general validity of the operator 
product expansion, or of the QCD sum rules formalism. We will instead deal with 
the well-known factl that instantons do spoil the operator product expansion by 
introducing corrections that are power-suppressed by 9 or more inverse powers of 
the momentum. In view of the large power suppression, it is clear that corrections 
of this type will behave more like a step function, below which the perturbative 
methods will certainly be inapplicable. In other words, we expect a value of the ratio 
m T /AqcD below which instanton corrections are of order one, and above which they 
are essentially zero. We point out that this is the only non-perturbative correction 
that can be computed explicitly in terms of the QCD parameter A and the quark 
masses alone, with no need of further phenomenological inputs. Computations of 
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this type of corrections have been carried out in refs. H and || in the context of 
e + e~ — > hadrons. Although intermediate formulae and results are the same in the 
two papers, the conclusions are different, since the authors of ref. conclude that 
instanton corrections are suppressed by four powers of the momentum, while in ref. || 
the effects are suppressed by twelve powers. For our purposes we are interested in the 
kind of result given in ref. 0, since the corrections mentioned in ref. should be 
viewed rather as instanton corrections to the value of the matrix elements of dimension 
four operators, which is usually determined empirically, and is not calculable from 
first principles. 

In ref. ||, it is claimed that instanton corrections in e + e~ — > hadrons become 
large for annihilation energies between 1 and 2 GeV. If a similar result held in the 
case of the r, it would spoil the analysis of ref. [0]. 

In order to obtain results that can be applied to our case we must complete the 
calculations of refs. || and || in three respects. First of all, we need to include the 
corrections to the axial current correlators, which were not considered there. Secondly, 
we must consider non-diagonal currents for flavours of different masses, which also 
were not considered previously. Lastly, we must perform an analysis in a definite 
subtraction scheme, in order to be able to relate the value of A that we intend to use 
in this context with the value that is extracted from high-energy experiments. In fact, 
instanton effects are suppressed by a factor of order exp(— 2-k/ols), so that a scheme 
redefinition, which changes the inverse of as by the addition of a constant, affects 
directly the prefactor of the instanton correction (at the time when refs. || and || 
were written, there was no sound agreement on the allowed range of Aqcd ; so that 
a complete answer correctly including the scheme dependence would not have been 
very useful). Our result also differs quantitatively from ref. [f|, and we will comment 
upon the differences in due time. 

Our paper is organized as follows: in section 2 we give a description of the calcu- 
lation of the process in question. In section 3 we review the formula for the instanton 
density in QCD, and specify its form in the MS scheme. In section 4 we apply our 
results to the r hadronic width. In section 5 we give our conclusions. 
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2. Instanton Corrections to the 

Current- Current Two-Point Functions 

We consider the correlator of two currents in the instanton background 
rV(x, y; Q ± ) = (j; d (x)J^(y)) A{ n ±) = - Tr (v,S d ± (x, y; Q ± )T v Sl{y, x\ fl±)) , (2.1) 
where we have 

j; d = jf = tp d T^ u ; (2.2) 

T M stands for 7^ for vector currents, and for 757^ for axial current; fl± denotes the 
instanton (anti-instanton) global coordinates Q± = (z,p,R), where z is the position, 
p is the size, and R stands for the colour orientation. S± denotes the fermion propa- 
gator in the instanton (anti-instanton) background, whereas So is the free fermionic 
propagator. The superscript u, d specifies the mass of the propagator. When no 
superscript is given, it is massless. 

We consider here the general case of non-diagonal flavour currents. The two 
flavours involved, which will be simply called u and d, have different masses m u and 
rrid- The propagators in the instanton background have the small mass expansion 



m 



+ m u{d) J dz 4 S ± {x,z;n±)S±{z,y;tt±) + 0{m 2 ), (2.3) 

where terms with odd powers of the mass commute with 75, while the terms with 
even powers anticommute. 

From eqs. (|2 . 1| ) and (|2.3|) we can immediately obtain the correlator in the limit of 
small masses 

u%, = -Tr(-y li S (x,y) 7v S (y,x)) + A^ V + CB, U (2.4) 
= -Tr(^S (x,y)j„S (y,x)) + A^-CB^, (2.5) 
where V stands for vector-vector and A stands for axial-axial correlators, and 
C - (— + md 



2 \m d m u , 

Aya, = Tr(-fyS (x,y)'y u So(y,x)) -Tr[y lx S±(x,y,Q±)'Y u S±(y,x;Q±)] (2.6) 



= 2Tr 



j^ (x)^l(y)^ v J d 4 zS±(x,z;tt±)S±(z,y; Q±) 



(2.7) 
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Expressions for the coefficients A^ v and B^ v have been obtained in ref . f| . Their result 
agreed with that of ref. || . We do however find a different colour factor normalization 
for the result, our normalization being 2/3 of theirs. We get (after colour averaging) 



A/wix, y, p) — -^iS^ayfi 



A 4 

P 



P\h x h y ) 2 (2A « A/3 _ g af3 A 2j 



+X-,h x h y (h y (A a y? + AV) - h x (A a x^ + A?x 



odd terms (2.8) 



1 p2 

Bf, v {x, y, p) = - — {Kh y ) 2 — [(p 2 + x ■ y)g^ + (y^x v - x^y u )\ + odd terms, (2.9) 
where 

h x = l/(x 2 + p 2 ), h y = l/(y 2 + p 2 ) 
A = x - y, 

"Odd terms" refers to terms that change sign when going from an instanton to an 
anti-instanton, and thereby vanish when summing over the two contributions. The 
instanton location has been chosen in x± = in the above formulae. 

The expression for the current- current correlator including instanton corrections 
in the dilute-gas approximation is 

rv(A) = n%(A) + / d A zd P D(p) £ (rv(*, y , z , P , ±) - nJ„(A)) , (2.11) 

where D(p) is the instanton density, which will be specified in the following section. 
Using eqs. (|2.4j) and (|2.5|) we obtain 

IV(A) = U%(A) + I d 4 zdpD(p)2 {A^{x -z,y-z,p) + CB^x - z,y - z,p)) 

(2.12) 

for vector currents, and the same form with C — > — C for axial currents. We have 

^ } ~ (2tt 2 ) 2 A8 ' 1 ' 

This formula includes the factor 3 from the colour trace. Defining 

a^(A,p) = J d 4 zA fll/ (x - z,y - z, p) (2.14) 
b^ v (A,p) = J d A zB^ u (x - z,y - z, p), (2.15) 



-5- 



the z integration gives 



2ir 2 



2tt 2 



d 2 
dA^dA 

d 2 



-G(A\p)+2G'(A\p)g : 



<9A< 



G{A\p) + 2G\A\p) 



(2.16) 
(2.17) 



where G and G' are functions of A 2 and p, and they are given by the expressions 



G\A\p) = 



dG(A 2 ,p) p 2 



dA 2 



A 4 



2p 2 1 , $ " 1 



"7 lo g 



A 2 £ °£+l 



(2.18) 



with £ = y 1 + 4p 2 /A 2 . The above expression is a non-analytic function of p for small 
p, while it is analytic in p' 1 for p — > oo: 



limG'(A 2 ,p) 



lim G'(A 2 ,p) 



A 4 



i 2/?2 1 ^ 
~ ~ A 2 " A 2 



1 1 

+ 



A 2 A 4 
+ 



6A 2 30p 2 140p 4 630p 6 2772p 8 



(2.19) 
(2.20) 



We need the integral of G against the instanton density, which is given by 

1 



D(p) = H 



log 



p 2 A 2 



P 



M 



(2.21) 



where M = 6 + n f /3. The value of H and of the power of the logarithmic term 
c will be discussed in more detail in the next section. We notice that the small-p 
region never poses a convergence problem. On the other hand, the large-p region is 
divergent. However, because of the analiticity of the integrand in that region, the 
divergence will be limited to a finite number of terms in the expansion of G. For 
definiteness, let us assume that rif = 3. Then D(p) will behave like p 7 for small 
p. Then we will have to subtract from and b^ u all the terms of their expansion 
up to the power p~ 8 . These subtracted terms will all have power-like dependence 
upon A 2 , with undefined, infrared divergent coefficients. We interpret these infrared 
divergent terms as the instanton contribution to the expectation value of operators in 
the operator product expansion of the two currents we are considering. These matrix 
elements are in general uncalculable, and they are usually estimated on the basis of 
some phenomenological considerations. The first term of this kind is in fact related 
to operators of dimension four. Since we want instead the true, direct effect of the 
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instanton upon our amplitude, we will subtract at the end these divergent terms. We 
first define 



&( A 2 , p, p ) = G'(A 2 , p) - 6(p - p )— £ 



i=o 



A 



2.; 



G(A 2 ,p,p o )=G(A 2 ,p)-0(p-po) 



logA 2 ^ + X:^f-) 



■2.1 



(2.22) 



^2.23) 



where the gj are the numerical coefficients in the expansion of eq. ( 2.20Q . The func- 
tions G, G' are then integrable in p in the whole range, even when multiplied by a 
power of p M with M < 2L — 1. Once we know that their integral is in fact conver- 
gent, we may regulate it in any way we like. For example, we may choose the analytic 
continuation method of taking —5 < M < —3, and then continuing to all the allowed 
values of M. We now compute the Mellin transforms. We get 



dpp M G'(£',p,p ) 



A M - 1 r(-M-4)r 2 



1 A A 2 Vo f " 2j+1 



M + 5 



sin 



Mvr 



A 



dpp M G(A 2 ,p,p ) 



, n M - 2j + 1 
2A M+1 



^2.24) 



M+l 



- log A 2 g 



r(-M-4)r 2 



M + 5\ /Mvr 



pr 1 



E 



M-2J+1 



(2.25) 



'M + l ^jM-2j + l 

Observe that the above formulae are now convergent for M < 2L — 1. The poles in M 
arising from the gamma functions cancel against those in the sums. From eq. ( |2.12|) , 
(gig ), (|2~T5| ), ( CT) and (gT?[) we get 



n A «/(A) = n° (A) - h 



log 



1 



C 



o 2 



dA 



'9^ 



d 2 



dA^dA 1 - 



A 2 A 2 
2A m+i 

M+l 



r(-M-4)r 2 



M + 5 



sm 



Mvr\ 1 



71 z 



(C-1)2^A A/ - 1 



D{p)dp— 

7T Z 



c 



a 2 



<9A : 



~9iiv 



d 2 



dAt*dA u 



logA 2 + (C-l)2^^ + P^(A 



(2.26) 



where P M „(A) is a polynomial in A. Observe that the logarithmic power in the 
instanton density is simply replaced by its value for p = p. Corrections to this 
replacement are logarithmically suppressed, and therefore are not included here. 
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The Fourier transform of the expression fl2.26| ) can also be performed for non- 
integer M, and then continued to the desired value of M, using the formula 



e iA ' p A J d 4 A 



P 



-J-4 



47T 



cos — 1 - 

2 V2 



r 



J + 3 



r(j + 4). 



(2.27) 



Observe that this formula vanishes when J is an even integer. In fact, in this case 
the result is a distribution concentrated at p 2 = 0. The term P jUV (A) therefore does 
not contribute to the Fourier transform, and we get the result 



3 I M+3 



(M + 3)r(|) r 3 
(M + i)r(^)" 



p 



-M-3 



+ 



A) 




p2 



+ (C - 1)^, 



3C)(^p 2 -^) + 3(l-C)p^ . (2.28 



The last term is explicitly infrared divergent for po — > oo- We interpret this term as the 
contribution of the instanton to the matrix elements of the dimension-four operators 
F 2 and mipip. There are no singularities in the remaining terms for positive values of 
M. 

It is easy to check that, in the particular case of flavour- diagonal vector currents, 
the l/p 4 term corresponds to the term obtained in refs. and 0. One can argue that 
the l/p 4 term cannot contribute to the discontinuity of the polarization operator, and 
therefore it does not affect the cross sections for hadron production. This formally 
correct argument fails however when radiative corrections are included. Similarly, we 
have neglected a number of contributions localized at p 2 = 0. Radiative corrections 
may delocalize these contributions, and therefore give corrections to the discontinuity, 
which behave like three or more inverse powers of p 2 . These terms cannot be neglected, 
and, being IR-divergent, cannot be computed by perturbative techniques. They are 
usually accounted for phenomenologically, when the value of the various condensates 
is extracted from data. 
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3. The Instanton Density 

The instanton density for SU(N) has been computed in ref. |7]]. In the absence of 
fermions it is given by the formula 



W 



PV 



4 exp[-a(l) - 2(N - 2)a(i)] f d 4 zdp /4tt 



x 



2iV 



cxp 



8tt 2 



9 2 (P) 



(JV-l)!(JV-2)! J p J 

where g is the strong coupling constant, and (see ref. ||) 

= 0.443307 

a i -j = 0.145873. 

The above formula is given in the Pauli-Villars regularization scheme, with 



(3.1) 



1 



57T 



9 2 ( P ) 



57T 



9' 



UN 



log(m p), 



(3.2) 
(3.3) 

(3.4) 



where m is the Pauli-Villars mass, and g without argument is the bare coupling. 
If we want to use the value of as measured in today's high-energy experiments, we 
should convert the above formula to the MS scheme. This change of scheme was first 
given in ref. M, and then corrected in ref. |J. We have also computed the change of 
scheme by just computing the vacuum polarization in the background gauge, in both 
the Pauli-Villars and the MS scheme. We find that the two schemes give the same 
result if the Pauli-Villars mass mo and the MS scale are related by the formula 



11 ml 1 
— log — - - 

3 /r 3 



0. 



(3.5) 



Expressing m as a function of \x and replacing it in formula ( |3.1| ) we obtain 



W 



MS 



4 exp 



-a(l)-2(iV-2)a(i)] rd'zdpfAn 2 



71 A 



2N 



exp 



8tt 2 



9 2 (P) 



(3.6) 



which agrees with ref. 
extra factor 



(JV-l)!(iV-2)! J pr \g* 

in the SU(2) case. Including the fermions, we obtain the 
2 



p ni Iiirai exp 



-m log(m p) + 2n f a (| 



(3.7) 



For the fermion contribution to the vacuum polarization, we find that the Pauli- 
Villars and MS results agree if p = m, so that in this case it is enough to replace 
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mo with p. We will also need to express the running mass in terms of invariant mass 
parameters 

/ 1 AM 33 ~ 2n f 
m[p) = m (log-r ) (3.8) 

according to the definition of ref . fij . Observe that consistency of the order at which 
we are computing requires that one uses the full two-loop expression for g in the 
exponent (this was not included in ref. ||) 



2tt 



Anbo log 



6(153 -19n f ) log log ^3 
+ (33-2n f ) 2 log^ 



(3.9) 



while in the prefactor a leading-order formula is accurate enough. Gathering all the 
factors, and setting p = 1/p, we get 



D(p) 



H 



H 



log 



p 2 A 2 



»£ 2 



(3.10) 



) A u+ 3 _ eX p [-a(l) + i + (2th - 2)a (i 



7T 



2 33-2n f 



/33-2r^ 6 
V 12~ 
45 - 5n f 



(3.11) 
(3.12) 



33 - 2m 

Our final formula for the instanton contribution to the vacuum polarization is then 



33 + nf 



u, u (p 2 ) = u°(p 2 ) + K [n 



rrij 

T 



p" 

A 2 



log 



V 

A 2 



X 



G (l0 + y) " l] (pV " PmP«) + (C - 1) (l0 + y) PmP-} (3-13) 



with 



— exp [-a(l) + | + {2n f -2)a (§)] 



7T 



33 - 2n { 
12 



2 33-2n f 



X 



c 



( 9 + f)r(|)r3(| + f ) 



7+f )r(6 + f 



±(m u /m d + m d /m u )/2 



(3.14) 



(3.15) 



where in the last expression the + sign is appropriate for vector-vector, and the — 
sign for axial-axial correlators. The Born term (which fixes our normalization) is 
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given by 

nj y = -^(PuPu ~ P 2 9nv) logp 2 (3.16) 

both for the axial and the vector contribution. 

All results quoted so far were obtained in the Euclidean metric. The corresponding 
Minkowski space formulae in the time-like region are obtained by analytic continua- 
tion in p 2 . Observe that in the case of the electromagnetic current, which has C = 1 
the sign of the instanton correction is opposite to that of the leading term, contrary 
to the result of ref. 0. 



4. Final Results 

It is now straightforward to obtain the instanton contribution to the r hadronic 
decays. According to ref. 0, using the same notation, the ratio of the hadronic to 
the leptonic width R T is given by 



R T = 6m f dz(l - zf [(1 + 2z)Tl%lv(P 2 ) + ngy(p 2 )] , (4-1) 

J\z\=l 

where z = p 2 /M 2 , and 

n£. v = Ti { Hv{P 2 WP V ~ pV) + rftlvtfW. (4.2) 
The suffix A + V indicates the sum of the axial and vector contributions. We get 

^r inst / A \ 9 m u rh d m s f , , „ fi / -zM 2 \^ r , 

(4.3) 

Observe now that the contour integral would give zero if the logarithmic term was not 
present. Since the power of the logarithm is very near 1, we expect that the integral 
will not depend much upon the ratio A/M T . In fact we find the numerical result 

R™ st ( 3. 64 A\ 9 m u m d m s 



i?o V Mr J Ml 



(4.4) 



where the coefficient 3.64 corresponds to A = 400 MeV; it varies from 3.67 to 3.62 if 
A is pushed to the extreme values of 0.1 and 1 GeV. Choosing with ref. [0 



rh u = 8.7 MeV, rh d = 15 MeV and rh s = 270 MeV 



(4.5) 



with M T = 1.784 GeV we get the result 



R r St -( 0-96 >< A \* (46) 
#° V 1.784 GeV/ 1 ' 1 

The value of A to be used in this context is A3, which is somewhat larger than the 
corresponding values of A5 that are usually quoted. For example, a recent review 
quotes the range 150 < A 5 < 330 MeV, which corresponds roughly to 280 < A 3 < 
510 MeV. We see that even in the most pessimistic case of the largest allowed value 
for A3 = 510 MeV the instanton correction would turn out to be absolutely negligible. 



Conclusions 



We have completed a calculation of the one-instanton contribution to the r hadronic 
width. We found that the corrections are actually negligible. Part of the smallness 
of the result is due to the fact that the one-instanton correction is proportional to 
the chiral suppression factor m u mdm s /M^. It is interesting to note that if that factor 
were not there, the instanton correction could be of order 1, and would thereby in- 
validate the conclusions of ref . [T] . At this point one may wonder whether corrections 
due to instanton anti-instanton pairs, which should be suppressed by 18 powers of 
the ratio A/M T , but do not carry any chiral suppression, may have coefficients of 
comparable size, that is to say, if they become of order 1 for M T 4A 3 . The com- 



putation of ref. [ 1 1] seems to indicate large corrections, although one may doubt the 
reliability of the method used there to perform a computation beyond the dilute gas 
approximation. 

Our result for the instanton correction to the vacuum polarization, for m — 3, can 
be summarized as follows 



n„„(p 2 



4tt 2 



(PvPv -P 2 9 nv) log P 2 + 



m u m d m s 
p 3 



/5.1701 A^ 


9 

) 


' p 2 ' 


V P / 







_, 10 

9 



X 



To) (P 2 9^u - P^Pu) + ts(C - l)ptf>„ 



with C = ±(m u / ma + rrid / m u ) . For electromagnetic currents (C 
the logarithmic exponent with 1, we get 



(5.1) 

1), approximating 
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with Q = 0.591 GeV for A 3 = 300 MeV and Q = 0.867 GeV for A 3 = 500 MeV. 
The considerable numerical differences with respect to ref. |J have several origins. 
First of all, there was an error in the normalization factor for the instanton density!!, 
which was subsequently corrected by the author himself. The authors of ref. use 
a leading-order expression for as (instead of a next-to-leading one) in the exponent 
of the instanton density, which leads to an overestimate of the effect. The other 
differences are due to the use of the MS instead of the Pauli-Villars scheme in the 
instanton density, the different definition of the invariant quark mass used in the 
present work, and the different colour factor we found. These last three differences 
have a minor impact on the final result. 
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